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SOME BI-MATRIX MODELS FOR BI-FREE LIMIT DISTRIBUTIONS 
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Abstract. In this paper, an analogue of matrix models from free probability is developed in the bi-free 
setting. A bi-matrix model is not simply a pair of matrix models, but a pair of matrix models where one 
element in the pair acts by left-multiplication on matrices and the other element acts via a ‘twisted’-right 
action. The asymptotic distributions of bi-matrix models of Gaussian random variables tend to bi-free central 
limit distributions with certain covariance matrices. Furthermore, many classical random matrix results 
immediately generalize to the bi-free setting. For example, bi-matrix models of left and right creation 
and annihilation operators on a Fock space have joint distributions equal to left and right creation and 
annihilation operators on a Fock space and are bi-freely independent from the left and right action of scalar 
matrices. Similar results hold for bi-matrix models of g-deformed left and right creation and annihilation 
operators provided asymptotic limits are considered. Finally bi-matrix models with asymptotic limits equal 
to Boolean independent central limit distributions and monotonically independent central limit distributions 
are constructed. 


1. Introduction 

Since the notion of bi-free pairs of faces was introduced by Voiculescu in m, the theory has quickly 
developed by generalizing many ideas and results from free probability. For example, |19| determined the 
bi-free central limit distributions. On the combinatorial side, Mastnak and Nica in [8] introduced a collection 
of partitions for bi-free pairs of faces that was postulated to be analogous to the role non-crossing partitions 
play in free probability. In [5], the postulate of Mastnak and Nica was confirmed to be correct. Subsequently 
[2] generalized such notions to the operator-valued setting where things get ‘interesting’. 

Following the development of free probability, Voiculescu in [20] constructed a bi-free partial i?-transform 
as an analogue of his i?-transform from m, which plays an important role in free probability. Furthermore, 
Voiculescu in |21j developed a bi-free partial S'-transform and combinatorial proofs were later found in 
m and [14] respectively. In addition, the notion of bi-free infinitely (additive) divisible distributions was 
developed in [5] thereby discovering the bi-free Poisson distributions. 

One important result in free probability is the connection between free probability theory and random 
matrix theory. Indeed Voiculescu in (TS] connected these two theories by demonstrating that the distributions 
of certain matrices tend to the free central limit distributions, namely semicircular distributions, as the size 
of the matrices increase. Furthermore, results pertaining to asymptotic freeness were demonstrated in [18] 
and similar results for other matrix models were developed in miiiiii]. The matrix models of HHIKIllTH] 
have since been important tools in free probability. 

The goal of this paper is to introduce a bi-free analogue of matrix models. Such bi-matrix models have 
been elusive and the view taken in this paper is that a bi-matrix model is not simply a pair of matrices, 
but a pair of matrices with specific actions on matrices. This approach of requiring actions on matrices may 
be slightly unsatisfactory, but the notion of random matrices acting on matrices has precedence (e.g. i) 
and work in m Section 6] provides evidence (beyond that presented here) for why this view is possibly the 
correct one. Including this introduction, this paper contains seven sections which are summarized as follows. 

The necessary background on bi-free probability is recalled in Section]^ Said background includes bi-non¬ 
crossing partitions, the bi-non-crossing Mobius function, (£, r)-cumulants, the universal moment polynomials 
for bi-freeness, and the abstract structures required for bi-freeness with amalgamation. Complete expositions 
of these results may found in W- 

The structures and notion of a bi-matrix model are introduced in Section [31 One should think of a bi¬ 
matrix model as a pair of matrices of operators that act on matrices of vectors (or operators) where the 
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left matrix in the pair acts via left matrix multiplication whereas the right matrix in the pair acts via right 
matrix multiplication with a ‘twist’. In addition to developing such notions, an essential lemma is presented 
which enables one to compute the joint distributions of bi-matrix models. 

Bi-matrix models with commutative entries are examined in Section 3] Due to the commutation of the 
entries, bi-matrix models simplify slightly in this setting. Many results analogous to those in free probability 
are obtained. In particular, bi-matrix models with asymptotic distributions equal to certain bi-free central 
limit distributions from m are demonstrated. However, due to commutative, only bi-free central limit 
distributions where the left and right operators commute in distribution may be obtained. These restriction 
produced by commutativity have already been seen in the theory of bi-free pairs of faces in that the operator 
model in [8] using left and right creation operators on a Fock space was insufficient to describe all bi-free 
distributions and a more complicated model was required (see 0). 

However, commutativity is pleasant in bi-free probability as it enables the importation of results from 
free probability. In fact. Theorem 14.131 demonstrates that asymptotic bi-freeness for commutative bi-matrix 
models is not much more than asymptotic freeness of random matrices. Consequently most (if not all) results 
relating random matrix theory to free probability should have an analogue in the bi-free setting. Therefore 
we only exhibit specific results such as realizing bi-free Poisson distributions via random pairs of Wishart 
matrices and demonstrating the asymptotic bi-freeness from constant matrices of random pairs of Gaussian 
(or Haar) unitary matrices. 

Some inadequacies from Section |4] are rectified in Section[5l In particular, generalizing [121 Theorem 5.2], 
bi-matrix models where the left matrices consists of left creation and annihilation operators on a Fock space 
and the right matrices consists of right creation and annihilation operators have joint distributions are equal 
to the joint distribution of left and right creation and annihilation operators. Using these models, all bi-free 
central limit distributions may be obtained. Furthermore, such bi-matrix models are bi-free from the left 
and right actions of constant matrices. 

Bi-matrix models involving g-deformed left and right creation and annihilation operators are examined 
in Section El In particular, asymptotic bi-freeness from block matrices is obtain and the asymptotic limit 
of these bi-matrix models is again left and right creation and annihilation operators on a Fock space. The 
results of this section generalize those of m to the bi-free setting and provide an alternate approach to 
Theorem 14 .11 1 bv the q = 1 case. 

Finally, in Section [3 bi-matrix models with asymptotic limits equal to Boolean independent central 
limit distributions and monotonically independent central limit distributions are constructed. These results 
require the use of products of left and right matrices to represent a single random variable in order to obtain 
the desired limit distributions. 

2. Background on Bi-Free Pairs of Faces 

In this section, we will summarize some essential combinatorial aspects of bi-free probability from [3] and 
the structures for operator-valued bi-free probability from [2]. We will only present the results essential to 
this paper and refer the reader to complete summarizes in [2] and [13] respectively. 

Bi-Preeness for Pairs of Faces. For the remainder of this paper, a map x ■ {!> • ■ • i "r} {£,r} is used to 

designate whether the operator in a sequence of n operators should be a left operator (when x(fc) = £) 
or a right operator (when x(fc) = r). Similarly, a map e : {1,..., n} —>■ itT is used to determine which colour 
from a collection K each operator is coloured. 

The main difference between the combinatorial aspects of free and bi-free probability stem from handling 
the following permutation. 

Definition 2.1. Given y : {1,..., n} —>■ {£, r}, if 

X“^(W) = {*1 < ■•• < *p} and x“^(W) = {Wi > ■•• > 
the permutation on {1,..., n} is defined by s^{k) = ik- 

Note the sequence (s;^(l),..., s^xin)) corresponds to, instead of reading {1,..., n} in the traditional order, 
reading X~^{{£}) (that is, the left nodes) in increasing order followed by reading X~^({’’}) (that is, the right 
nodes) in decreasing order. 

Let 'P{n) denote the set of partitions on n elements. Given a partition tt G V{n) and x,y € {1,..., n}, 
we write x y whenever x and y are in the same block of tt (a set in tt) and x y otherwise. Given two 
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partitions tt, cr G 'P{n), tt is said to be a refinement of cr, denoted tt < tr, if every block of tt is contained in a 
single block of a. Refinement defines a partial ordering on V(n) turning T’{n) into a lattice. 

Definition 2.2. A partition tt G 'P{n) is said to be bi-non-crossing with respect to x if tfi® partition s~^ ■ tt 
(the partition formed by applying to the blocks of tt) is non-crossing. Equivalently tt is bi-non-crossing 
if TT is non-crossing on {1, ... ,n} when the set is ordered via Sj^(l) s^(2) s^(ji). The set of 

bi-non-crossing partitions with respect to y is denoted by BNC{x)- 

Example 2.3. If y : 6} —>■ {^,r} is such that X~^i.{^}) = {2,3,6} and = {1,4,5}, then 

(s;^(l), ..., S;^(6)) = (2, 3,6, 5,4, 1). If tt = {{1, 4}, {2, 5}, {3,6}}, then tt is crossing on {!,..., 6} yet is 
bi-non-crossing with respect to y. This may be seen via the following diagrams. 



One function on pairs of elements of BNC(x) required in this paper is the following. 

Definition 2.4. The bi-non-crossing Mobius function is the function 

piBNC ■■ U U BNC{x) X BNCix) ^ c 

defined such that p.BNci'^, <^) = 0 unless tt is a refinement of tr, and otherwise defined recursively via the 
formulae 

X] h‘BNc{j,a)= ^ AiBArc(7r,T) = 

reBNCix) reBNCix) 

7r<r<<7 7r<r<(T 

Due to the similarity in lattice structures, the bi-non-crossing Mobius function is related to the non¬ 
crossing Mobius function ^atc by the formula 

AtSArc(7r, cr) = UNcis^^ ’ tt, • a). 

This implies that fiBNC inherits many ‘multiplicative’ properties that p,NC has. For more details, see 
[21 Section 3]. 

Definition 2.5. Let (M, (p) be a non-commutative probability space: that is, let M be a unital algebra and 
M —>■ C be unital and linear. A pair of faces in A is a. pair (C, D) of unital subalgebras of A. We will 
call C the left face and D the right face. 

Given Zi,..., Zn G A and tt G V{n) with blocks 14 = {fci.a, < • • • < for x G {1,... ,p}, define 

p 

PniZl, . . . , Zn) '■= 

31—1 

Given a map y : {1,..., n} —>•{£, r} and tt G BNC{x), define 

CTtt , ■ • ■ , Zyi) . — ^ ^ Pa , ■ • ■ , Zji)pjQjXfQ (iL, tt) . (1) 

aeBNCix) 

(T<7r 

The {£, r)-cumulant of Zi,..., Zn corresponding to y is k^{Zi, ..., Z„) := {Zi, ..., Zn) where is the 
full partition {{1,..., n}}. Each (£, r)-cumulant should always be viewed as a function where only elements 
of left faces may be inserted into the lA^ entry when y(fc) = t and only elements of right faces may be inserted 
into the fc**' entry when y(A:) = r. Furthermore, if tt G BNC{x) has blocks 14 = {ki^x < ■ ■ ■ < for 

a; G {1,..., p}, then one can verify 

p 

CTtz^Zi , . . . , Zn) ^'^1 Va; ■ 5 .3: )■ 

x—1 


( 1 if TT = a 
\ 0 otherwise 
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Furthermore, the reversion formula 

'^{Zi---Zn)= ^ K^{Zi,...,Zn) ( 2 ) 

n&BNC(x) 

holds thereby giving each moment as a sum of products of {£, r)-cumulants. 

The following is the main result from [3]. For e : {1,..., n} —?> iF, note e defines an element of 7^(n) whose 
blocks are {e~^{{k})}k^K■ 

Theorem 2.6 ([3] Theorem 4.3.1]). Let {A, if) be a non-eommutative probability space and let {{Ck^ Dk)}k&K 
be a family of pairs of faces from A. Then {(Cfc, Dk)}keK are bi-free with respect to ip if and only if for all 
X : {1,..., n} r}, e : {1,..., n} -)> K, and 

*/x(w) = ^ 

ifxi'ni)=r ’ 


Zm € 




(m) 


have 


<^(Zi • • • Z„) = ^ 

Tr&BNCix) 


^ MBArc(7r,cr) 
aGBNC(x) 

7r<(7<e 


^TT (■^1 ^n') • 


(3) 


Furthermore, it suffiees to verify equation m for Zk in generating subsets of Cf^(^k) and D^(^k) and bi-freeness 
of {{Ck, Dk)}k€K is equivalent to k^[Zi, ..., Zn) = 0 whenever x, £, and Zk are as above and e is not 
eonstant. 


Important to this paper are the bi-free central limit distributions of m Section 7]. 

Definition 2.7. Let (.4, (p) be a non-commutative probability space and let I and J be disjoint index 
sets. A two-faced family {{Zi}i^i, {Zj}j^j) is said to be a (centred) bi-free eentral limit distribution if all 
{£, r)-cumulants of order 1 and of order at least 3 are zero. 

Note [m Theorem 7.4] completely describes all possible (centred) bi-free central limit distribution via the 
(complex) numbers ip{ZkiZk^) = K^(Zkj^, for ki,k 2 € I U J, and concretely represents such two-faced 
families as combinations of left and right creation and annihilation operators on Fock spaces. The matrix C 
indexed by / U J obtained with the (fci, ^ 2 ) entry equal to ip(ZkiZk^) is called the eovariance matrix of the 
family ({ZJ^g/, 

Structures for Operator-Valued Bi-Fheeness. With [3] things just got stranger as there are two copies 
of the amalgamation algebra and the natural structures for operator-valued bi-free probability appear, on 
the surface, very different from from those for operator-valued free probability. We describe these structures 
here. For the remainder of this section B will denote a unital algebra. One may replace B with the N x N 
matrix algebra A4jv(C) as this is all that is required in other sections of this paper. 

Definition 2.8. A B-non-commutative probability space is a pair (A, $) where A is a unital algebra con¬ 
taining B (with lx = Is) and $ : A —>■ B is a unital linear map such that 

$(&lZ&2) = bi^{Z)b2 

for all &i, &2 G B and Z G A. 

Definition 2.9. A B-B-non-eommutative probability spaee is a triple {A, Ej^, e) where A is a unital algebra, 
£ •. B ® B°'k‘ —>• A is a unital homomorphism such that £\b®1b ^nd ejiB^sop are injective, and Ej( ■. A^ B 
is a linear map such that 

EA{£{bi®b 2 )Z) = biEj\^{Z)b 2 and Ax(^e(& 0 Is)) = Ax(2'e(lB (g) &)) 

for all bi,b 2 ,b G B and Z G A. To simplify notation, Lb and Rb are used in place of £{b 01b) and e(ls 0 &) 
respectively. 

The unital subalgebras of A defined by 

Ai := {Z G A 1 ZRb = RbZ for all b G B} and 
Ar '■= {Z G A 1 ZLb = LbZ for all b G B} 
are called the left and right algebras of A respectively. 
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Given a i3-_B-non-commutative probability space (^, s), notice {Ai, Ej\) is always a i?-non-commutative 

probability space with e{B ( 8 ) 1b) as the copy of B and {Ar, EX) is a i3°P-non-commutative probability space 
with £(1 b < 8 ) i?°P) as the copy of B°'^. Furthermore {A, Ej\^, e) is simply a non-commutative probability space 
whenever B = <C. 

The reason _B-i?-non-commutative probability spaces are the correct setting of operator-valued bi-free 
probability can be seen through the following. 

Definition 2.10. A B-B-bimodule with a specified B-vector state is a triple (A, X,px) where A is a direct 
sum of B-H-bimodules 

X = B®X, 

and px : X ^ B is the linear map 

Pxib®r]) = b. 

Let £(A’) denote the set of linear operators on X. For each b G B dehne the operators Lf,, Rf, G E(X) by 

Eb(v) = b - r] and Rb{v) = V ’ b for all rj G X . 

The unital subalgebras of C{X) defined by 

Ct{X) ■= {Z G C{X) I ZRb = RbZ for all b G B) and 

Cr{X) := {Z G C{X) I ZLb = LbZ for all b G B} 

are called the left and right algebras of L{X) respectively. 

Given a B-H-bimodule with a specified i?-vector state {X,X,p), the expectation of C{X) onto B is the 
linear map —>• B dehned by 

Ec(x){Z) = px{Z{Ib)) 

for all Z G C{X). 

Notice if £ : B G B°p —^ C.{X) is defined by s{bi 0 62 ) = Lb^Rb^i then {C[X), Ec(x)t^) is a concrete 
S-H-non-commutative probability space. Moreover [2l Theorem 3.2.4] demonstrated every abstract B-B- 
non-commutative probability space can be represented inside a concrete B-H-non-commutative probability 
space, at least for the purposes of computing distributions. 

3. General Structure for Bi-Matrix Models 

The general structures required for bi-matrix models are introduced in this section. These structures 
are motivated by the main result of m Section 6 ] which is stated below. The symbols 
denote the standard matrix units for MnX)- Either [a^j] or ® used to denote the 

N X N matrix with {i,j)^^ entry Oij depending on which better fits the current context. Furthermore Tr 
will denote the trace on A4Ar(C) defined by 

N 

i=l 

The following introduces the general structure for bi-matrix models for matrices with elements in C{X) 
where A is a pointed vector space. 

Construction 3.1. Let {X,X,^,px) be a pointed vector space; that is, A is a vector space over C with 
X = X and px ■ X ^ C is the linear map defined by 

Pa (AC © 77) = A. 

For N gN consider the A4Ar(C)-A4Ar(C)-bimodule Xn ■= M.n{X) where 




■ N 



■ N 

Ujj] ■ 

Put] = 

.k^l 

and 

{Vi,j] ■ — 

^ ^ ^k,j‘ni,k 
.k^l 


for all [oij] G MnX) und G Xm- Then X^ becomes an A4iv(C)-Af Ar(C)-bimodule with specified 

A4iv(C)-vector state via 


Xpf = 
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and the linear map pxm '■ —>■ A4n (C) defined by 


PX;,(h,j]) = [px{v^,j)]■ 


We call X]y the n{'C.)- bimodule associated with {X,px) and {£{Xn), Ec(^Xn)) the Mn{£^)-M.n{£^)- 

non-commutative probability space associated with {X,px)- 

Since A” is a pointed vector space, there is a canonical choice of linear functional (px ■ £{X) —> C defined 
by 

px{Z)=px{Z{0)- 

Therefore, the expectation Ec{Xn) ■ £{Xn) AfAr(C) is defined by 

Ec{XM)iZ) = px^iZlN,^) 


where In,^ is the diagonal matrix diag(^, ..., 

To consider bi-matrix models, two homomorphisms from MN{£iX)) into C{Xn) are required. 
L : Mn{C[X)) ^ C{Xn) by 




■ N 

Ti^k{pk,j) 

.k=l 


Define 


for all [r]ij\ S Xm and \Tij\ G A4 n{£{X)). Is is elementary to verify that L is a homomorphism, that 
L{Mn{£{X))) C £^{Xn), and that L{[aijlx]) = j]- On the other hand (or perhaps, on the other face), 

define R : 7Wjv(>C(A’)°p)°p ^ £{Xn) by 




■ N 

^ ^ Sk,jijli,k) 
.fc=l 


for all [pi^j] G Xm and G M.m{£{X)). Is is elementary to verify that i? is a homomorphism (hence the 
ops), that R{Mn{£{X)°'^)°'^) C £r{Xi\i), and that R{[aijlx]) = R[ai,j]- We use L and R instead of Lw 
and Rm as the size of the input matrices determine N. The image of L (respectively R) is called the left 
(respectively right) matrix algebras of £{X) and elements of this algebra are called left (respectively right) 
matrices of £{X). 


Remark 3.2. If G Mn{£[X)) are such that TijSk,m = Sk^mTij for all i,j,k,m, then it is 

elementary to verify that L([Tij])R{[Sij]) = R{[Sij])L{[Tij]). 

Remark 3.3. It is worth mentioning that for all [Zij] G A4 n{£{X)), 

Ri[Zij])lN,^ = R{[Zi,j])lN,^- 

This fact will be important later in the paper in order to ‘transition from right operators to left operators’; 
a theme that has been essential in bi-free probability (see [11I31IS])- 


Remark 3.4. Using the notation as in Construction 13 .1 1 if [Zij] G A4n{£{X)) then 

Ec^x^^^{L{[Z,,,]))=PxAmz^,J])INx)=PxA[Z^AO]) = [pxiz.AO)] = [px{Z.A]=Ecix.)mz^,j])). 

As similar computations hold for any product of left and right matrices of £{X), we see that E^(^Xm) applied 
to a product of left and right matrices of £{X) depends only on the sequence of left and right matrices, the 
entries in the left and right matrices, and the linear map px '■ E{X) —> C. In particular, for such products, 
Ec{Xn) is fbe expected expectation obtained by applying px to each entry of an element of M^n{£{X)). 


Remark 3.5. Given a non-commutative probability space (A, ip), there are many ways to view A C £{X) 
in a state-preserving way. In particular, we adopt the convention that X = A with px{A) = ip{A) (so 
= C © ker(y))) and A C £{X) via T{A) = TA for all A G X and T G A C £{X). In particular, this 
convention implies 

^Ar(T) - PxiniA)) = Px{T) = p{T). 

Under this convention and using the notation of Construction 13.li we see for e N that = A^Ar(v4) and 




■ N 

-k^l 


and R{[Zij])[Aij] 


■ N 

^ ^ Zk,j 

-k^l 
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for all [A-ij] G A4iv(A). This computation shows why one might call i? a ‘twisted’-right action as it 

looks like multiplying the matrices [Aij][Zij] in Mn{A) except that the opposite multiplication on A is 
used in each matrix entry. 

Using the same idea as in Remark l^l^ one sees that the joint distribution of L(A^jv(^)) and 
depends only on the sequence of left and right matrices, the entries in the left and right matrices, and the 
linear map ^ > C. 

It is important to emphasize one does not want the operators Zi^ in R([Zij]) to be elements of A acting 
‘on the right’ of A; that is, via the action S'(A) = AS. One reason for this, other than Theorem 13.61 is that 
L{[Tij]) and would commute by Remark 13.21 and too much commutativity in bi-free probability is 

not optimal. 

Previous evidence towards why Construction 13.11 is the correct mathematical construction to use is the 
following result which demonstrates how bi-freeness of pairs or faces extends to bi-freeness of matrices of the 
pairs of faces. 

Theorem 3.6 (specific case of |131 Theorem 6.3.1]). Let (.A, (/j) be a non-commutative probability space and 
let {{Ck, Dk)}keK be bi-free pairs of faces with respect to (p. For each N € N let {C{Apf), Ec{Xff)) be the 
canonical M.N(fC)-M.]si{C)-non-commutative probability space associated with A as in Remark \8.A Then 
{(L(Ad7v(C'fc)), i?(Ad7v(T*fc^)°^))}/ceAr bi-free with amalgamation over A4]s[{C) with respect to Ec{Xn)- 

A piece of the proof of Theorem 13.61 essential to this paper is the following computation. 


Lemma 3.7. Using the notation and conventions of Construction HOI let x ■ ^ and let 

Zk = L{[Zij.^k]) ifxik) = £ and Zk = R{[Zij.^k]) if x{k) = r. Then 

N 

(Zi ■ ■ ■ Z„)(I]Y^^) = ^ ,ji;l O ■ ■ ■ O C) E^[[ii , . . . , 1„), (ji, . . . , J„); Af) 


where 


5 ■ • ■ ) *n)) (jl) • ■ • ) Jn); At) — E,. 




(iV)---A,, 


Pn),3s^(n) 


(IV) e M^(C). 


Proof. By the linearity of all maps involved, it suffices to consider Zk = L{Zi^j^ Ei^j^(N)) when x(fc) = £ 
and Zk = R(Zi^^j^ ® Ei^,]k(AI)) when x(fc) = r. Note 

Zn{lN,{) = 0 Ei^j^(N) 

regardless of the value of x(n). 

If x(n — 1) = i, observe that 


Zn-lZn(lN,i) = ® Ei^_^j^_^(N)Ei^^j^(N) 

whereas, if x(n — 1) = r, observe that 

Zn-lZn(lN ,£,) = i ,i„_i (?)) <8 A’in (-^) j ^ (A^) . 

In particular, this pattern repeats where Eij^^jj^(N) is placed on the left-hand side of the product of matrix 
units if x(fc) = £ whereas (N) is placed on the right-hand side of the product of matrix units if x(fc) = r. 
Hence the result follows by the definition of s^. □ 


Remark 3.8. For Ey^((ii, ..., in), (ji,... ,jn)',N) to be non-zero and to be a diagonal entry (thereby con¬ 
tributing to the trace), it is required that jg (^k) = is (fe-i-i) foi' ^ (where n -I- 1 —^ 1). This implies 
A = = U,(.-Rfe)+i) for all k. 


4. The Commutative Case 

In this section, we will study bi-matrix models in the commutative case to obtain (commutative in distri¬ 
bution) bi-free central limit distributions. In particular, for this section we take A = (LooiVl, fi), E) for our 
non-commutative (well, commutative) probability space where 

E{f)= f fix) dpi. 

Jn 
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Although some results in this section may be obtained via the bi-matrix models with bosonic creation and 
annihilation in Section [S] (i.e. with g = 1), the ideas in this section are not only simpler but important for 
demonstrating that many random matrix models have direct analogues in the bi-free setting. In particular, 
this section can be summarized as, “there is a bi-free analogue of any random matrix model that can be 
generalized from random matrices with independent entries to pairs of random matrices where each pair is 
independent from each other pair and each pair has a certain covariance matrix associated to it.” 

We begin by recalling some definitions. 

Definition 4.1. A family Xi ,..., Xn of self-adjoint random variables in (Loo(f^, m), E) is a (centred) Gauss¬ 
ian family if there exists a non-singular positive n x n matrix C with real entries (called the covariance 
matrix) such that for all A: G N and all 1 < fi, ..., < n, 

E{X,, ■ ■ ■ X,^) =—=^=== [ dxi---dxn 

y (27r)" det(C) Jr" 

where x = (xi,..., x„) and (•, •) denotes the standard inner product on K.”. 

A family of complex random variables Ai,...,A„ in a ^-probability space (Loo(n,/r),A) is a complex 
Gaussian family if Re(Ai),..., Re(A„), Im(Ai),..., Im(Xji) is a Gaussian family. 


Of important use in this section is the following formula. 

Theorem 4.2 (Wick’s Formula; see [101 Theorem 22.3] for example). Let Xi,..., A„ be a Gaussian family 
in fj,), E) with covariance matrix C = [c^j]. For all k and 1 < fci ,... ,k„ < n, 

E{Xk,---Xkj= Y. n 

77 ^ 7^2 ( ti ) {x,y}€n 


where 7^2 (^) denotes all pair partitions on {1,..., n}. Furthermore 

E{XiXj) = Cij. 


Example 4.3. If 


Q _ Cx,X Cx,Y 
[ Cy,X Cy,Y 

is a non-singular positive matrix with real entries, then, by Definition 14.11 and Theorem 14.21 there exists a 
Gaussian family X, Y such that 


E{X^) = cx,x, E{XY) = cx,Y, E{YX) = cy,x, 


and E{Y'^) = cy,y- 


Example 4.4. If 


(J _ Cx,X Cx,Y 
I Cy,X Cy,Y 

is a non-singular positive matrix with real entries, then we claim there exists a complex Gaussian family 
X, Y such that 

E(X^) = 0, E(XY) = 0, E(YX) = 0, and E(Y^) = 0 


yet 

E (XX) = cx,x, E (XY) = cx.F, E (YX) = cy,x, and E (YY) = cy,y. 
Indeed consider the Gaussian family Zi,..., Z 4 with the covariance matrix 


i(C®L) = i 


Cx,x 

0 

CX,Y 

0 

0 

cx,x 

0 

CX,Y 

Cy,X 

0 

cy,y 

0 

0 

Cy,X 

0 

cy,y 


which is clearly a non-singular positive matrix with real entries. Therefore, if A = Zi+iZ 2 and Y = Z 3 -|-fZ 4 , 
then A, T is a complex Gaussian family that satisfies the above equations. Indeed 


E (a2) = E (Zl) + iE(ZiZ2) + iE{Z2Zi) - E (Z|) = 0, 

E (XX) = E (Zi^ Y iE(ZiZ2) — iE(Z2Z\) Y E (-^ 2 ) “ 
E(XY) = EiZiZs) Y iE(Z2Z3) Y iEiZ^Zi) - ^(^ 2 ^ 4 ) = 0, 


E (AT) = E{ZiZ 3) Y iE{Z2Z3) - iE{ZiZi) Y ^(^ 2 ^ 4 ) = cx,y, 
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and similar computations yield the other equalities. 

We now introduce Construction o into this setting. 

Definition 4.5. For N G N a,n N x N random pair of matrices on {Loo{^, fr), E) is a pair where 

X^ is a left matrix and X^ is a right matrix with entries from Loo(r2,/r) G fJ-))- 


Remark 4.6. Note that a random pair of matrices is not simply a pair of random matrices, but a pair of 
random matrices with a certain action on elements of Mn{Loc(S^, fEj). 


Definition 4.7. Let 


C = 


Cl,r 

C’P ^ C'p 'p 


be a non-singular, positive matrix with real entries. A self-adjoint C-Gaussian random pair of matrices is an 
N X N random pair of matrices (A^, A’’) on {Looi^, fJ^), E) with A^ = L{[XE]) and A*" = R{[Xlj\) where 

(1) AA G Loo(fi,M) for all k G {t,r} and i,j G {1,..., A}, 

(2) AA = Xj ^ for all k G {£, r} and i, j G {1,..., A}, and 

(3) {aC, Re(AA), Im(AA) | k G {i,r},i,j G {1,... ,N},i < j] is a Gaussian family such that 




Remark 4.8. Given a non-singular, positive matrix with real entries 


C = 


ce,e ce^r 

Cp^p 


one can always construct an A x A self-adjoint C-Gaussian random pair of matrices by Definition 14.11 and 
by taking direct sums of the covariance matrices from Examples 14.31 and 14.41 


Remark 4.9. Due to commutativity and positivity requirements, only certain bi-free central limit distri¬ 
butions {Zj}j^j) on {A,ip) will be limits distributions of random pairs of matrices. Indeed we 

will require that [A, <p) is a *-non-commutative probability space, is a positive linear functional, and the 
covariance matrix of ({Zi}i^i, {Zj}j^j) is positive, non-singular, and has real entries. In particular, each 
Zi and Zj must be a non-zero self-adjoint semicircular variable. Furthermore, as the covariance matrix 
is self-adjoint with real entries, K^{Zk-^, Z^^) = ^fci) for all ki,k 2 G / U J. Since a joint moment 

of {{Zi}i^j,{Zj]j^j) can be computed as sums of products of second-order (.^, r)-cumulants by the defini¬ 
tions and results of Section [2l one can verify that {{Zi\i^i,{Zj'\j^j) commute in distribution. If C is the 
covariance of ({Zi}i^i,{Zj\j^j), we call {{Zi}i^i,{Zj}j^j) a C-bi-free central limit distribution. 


Theorem 4.10. Let 


C = 


C£^r 

Cp^£ Cp^p 


be a non-singular, positive matrix with real entries and let {Si, Sr) be a C-bi-free central limit distribution 
with respect to if. For each N G N let (A^(A), A’’(A)) be an N x N self-adjoint C-Gaussian random 
pair of matrices. Then the joint distribution of (A^(A), A'’(A)) with respect ;^Tr o Ejr(^;^;^) tends to the 
joint distribution of {Sg,Sr) with respect to if as N tends to infinity; that is, for every n G N and every 
X n} {l,r}. 


= If . 
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Proof. The proof presented here is motivated by the proofs in [TOl Lecture 22]. Notice, by Lemma 13.71 
Remark 13.81 and Theorem 14.21 that 


1 

N 


Tr(^Ecix) 


1 


N 


ATf + l 

1 

fVt +1 

1 

Nt+i 


E rp ( vXll) -^X( 2 ) ... V 

-^ 1 ^ 21,2 1 ^ 22,2 1 ‘ ‘ I ,-l 

\ sv(sv (1) + 1) ®v(sv (2) + l) 


2 l,...,2n = l 

N 

i; i: n 

2 l,...,2n —1 7TGV2(n) {x,y}G'7T 

N 


x(") 




x{y) 


’xUx^U) + i) ^xUx^(^/) + ^) 


E s n ^ N'l..-.,.,*. 


il,-..lin —1 Tre'P 2 (n) {Xjyy^TT 


X 


xiv) 

u.-i, 




(4) 

(where the last line follows by replacing ik with Notice equation (j4]) is zero unless n = 2m, in which 

case it equals 

N 

7rGp2(2m) ji,.” J 2 m = l {x,y}G'rT 
N 

EE n ,jy+l ^jy ,jx + l Cx(Sx (x)),x{Sx iv)) 

7re2’2(2m) ii ,...,t2m = l {sx(a:),Sx(?/)}6’^ 

N 

EE n JH+l^iHUx+lCx(sx(a;)).x(sx(y))' 


1 

JSfi+ri 


1 

_/Vi+" 

1 


7\ri+" 

7re'P2(2m) ti,...,i 2 m = l {x.ytgSx^'"^ 

The computations on [101 page 365] demonstrate for a fixed tt G V2{‘2m) that 

1 ^ 

iVi+'" E n '^ix.i^+i<^i,;jx+i = 0 

tl, - U2m = l {x,y}eTr 

unless TT is non-crossing in which case the limit is 1. Consequently, for a fixed tt G 7^2 (2m), 

N 


1 


lim 

A-)-oo 


E n '^*x,jH+l'^iy.ix + lCx(. 


Sx(^))A(sx(y)) 


= 0 


(5) 


^iv5'^2m —1 {ai,y}Gs 5 j; ^-Tr 


unless s“^ • tt is non-crossing pair partition, which means tt G BNC(x) is a pair partition. For a pair 
partition tt G BNC{x), the limit in equation ([5]) is Hjx y}e-K'^xix),xiy)- Therefore, if BNC^ix) denotes the 
bi-non-crossing pair partitions corresponding to y, then 

a! 1 "L (v>=<->(iV)--.Y* 


~ E n ^x(x),x(v) 

'K^BNC2{x) {37,y}€'7r 

= ^ K,r(<S'j(.(l), S'x( 2)5 ■ • ■ ,'S'x(n)) 

tt^BNC2{x) 

~ P (‘^x(i) ■ ■ ■ '^x(")) • 


□ 


Using the above ideas, asymptotic bi-freeness of random pairs of matrices is easily obtained via observing 
the correct ‘colouring’. 


Theorem 4.11. Fix an index set K. For each k G K let Ck be a2x2 non-singular, positive matrix with real 
entries and let {{Six^ Srx)}k^K be a collection of bi-free two-faced pairs with respect to if where (SixjSrx) 
is a Ck-bi-free central limit distribution. For each iV G N and k G K let be an N x N 

self-adjoint Ck-Gaussian random pair of matrices such that entries are independent for different k G K. 

For every n G N, every x ■ {Ij • ■ •) ''t-} {^j ^}j every e : {1,..., n} —^ RT, 

W)) = If (5,(1),(1) • • • ^x(n),dn)) ■ 
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In particular, are asymptotically bi-free with respect o Ec(Xf^). 

Proof. By repeating the ideas of Theorem 14.101 we obtain 
1 

N 


pTr (^Ecix) 


1 


N 


E E n 




jVt + l ^ ^ 11 ^E-\.)’E-\y) + lis-\y)’E-\.) + l^(^)^^(vrxix),x(v)- 

Tr€'P2(n) ii,... Jn = l {x,y}GTr 


The result now follows as in Theorem 14.101 as the only tt that contribute asymptotically to the sum are 
TT S BNC 2 {x) that may be coloured correctly. □ 


Remark 4.12. A more eeneral result than Theorem 14.111 is easv to obtain. Indeed let fiZ.Tcr. 

be a (centred) bi-free central limit distribution and let C = {ck,m)k,meiuj be the (/U J) x (/U J) covariance 

matrix. 

If C is non-singular, positive matrix with real entries, then similar arguments to those given in Remark 
l4.8l show that for each N gN there exists {X^j{N) \ i,jG{l,..., N}, k G lUj} C (Loo(fl, ti), E) such that 

• XX(^N) = Xj(j(iV) for all fc G / U J and i, j G {1,..., N}, and 

• I i, j G {I,N},i < j,k G I Li is a Gaussian family such 
that 

E[xl]iN)X^UN)) = 

For each N gN and A: G / U J, let X^{N) = Z([AA]) where Z = L ii k G I and Z = R\ik G J ■ Similar 
arguments to those above show that the asymptotic joint distribution of {{X^{N)}i^j,{X^{N)}j^j) with 
respect to ;^Tr o Ec{Xn) as —>■ oo is equal to the joint distribution of {{Zifi^i, {Zj}j^j). Theorem 14.111 
then follows by using the correct covariance matrix. 


For those familiar with the development of bi-free probability, it should not be a surprise that results 
from free probability can be generalized to bi-free probability provided all left operators commute with all 
right operators. Indeed the following result can be interpreted as “in the commutative world, asymptotic 
bi-freeness is pretty much asymptotic freeness.” 

Theorem 4.13. Let K be a fixed set and consider {C{X),ipx) for some pointed vector space X. For each 
N G N let {Xk{N)}k^K be left N x N matrices and let {YklN)}k£K be right N x N matrices of C{X). 
Suppose that the joint distributions of {Xk{N)}k£K U {Yk{N)}k^K converge with respect to ;^Tr o E^^n)- 
Furthermore, suppose 

(1) Xk{N)Ym{N) = Yjn{N)Xk{N) for all k,mG K and for all N, and 

(2) XklN)lN = Yk{N)lN for all k G K and N. 

If K = K 1 UK 2 then {Xk{N)}k^Ki is asymptotically free from {Xk{N)}k^K 2 with respect to ^XvoE c,{Xm) */ 
and only if {{Xk{N)}k^Ki,{Yk{N)}keKi) is asymptotically bi-free from {{Xk{N)}k^K 2 , {YkiN)}keK 2 ) with 
respect to jfYr o Ec(Xn). 

Proof. The proof of this result is obtained by demonstrating that the universal bi-free moment polynomials 
(i.e. equation dS])) hold asymptotically for {{XkiN)}keKiAYk{N)}keKi) with {{Xk{N)}keK 2 ,{YkiN)}keK 2 ) 
if and only if the universal free polynomials (those in equation ([3|) restricted to y with xi'OT-) = ^ for all m) 
hold asymptotically for {Xk{N)}k^Ki with {Xk{N)}k£K 2 - carbon copy of the proof of |H Theorem 10.2.1] 
with lim^v-^oo inserted in multiple places yields the result. We omit further details as similar arguments will 
be used later in the paper (see Example 14.151 and the proof of Theorem 15.31) . □ 


Remark 4.14. Note the condition ^^{Xk{N)}keKi is asymptotically free from {Xk{N)}k£K 2 with respect 
to ;^Tr o Ec(Xn)’' in Theorem 14.131 is precisely saying that {Xk{N)}k^Ki is asymptotically free from 
{Xk{N)}k^K 2 as usual matrices of operators. As such, since random matrices have commutative entries, 
many results from random matrices immediate have random pairs of matrices analogues. Furthermore, if a 
result can be proved for random matrices grouped into pairs of two, one immediately has a random pair of 
matrices result by viewing the first element in each pair as a left matrix and the second element in the pair 
as a right matrix. 
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Example 4.15. As an example application of how Theorem HUS] and its proof work, we will demonstrate for 
every bi-free Poisson distribution /i from 0 Example 3.11b] a random pair of matrices with limit distribution 
Fix a number A G (0,1) and for each N G N choose Mjq G N such that limAr-s-oo = A. For each 
G N let Y{N) be an Af x Mm random matrix whose entries are independent Gaussian random variables 
with mean 0 and variance 1. If 

X{N) = -^Y{NYY{N) 

(a Wishart matrix) then the Marchenko-Pastur law implies the limit distribution of X{N) with respect to 
■YTt o Ec(Xpf) is equal to the free Poisson law fip with rate A and jump size 1. Recall the free cumulants 
of p-p is A for all n G N so the free cumulant of X{N) with ;^Tr o tends to A as A^ tends to 

inhnity. 

For each {a,/S) G [5J Example 3.11b] defined the bi-free Poisson distribution with rate A and jump 
size {a, 13), denoted ptbP,a,/3, to be the limit distribution of 

// \\ \ \ 

where fflffl denotes the bi-free additive convolution. Furthermore [5] demonstrates that if y : {1,..., n} —>■ 
{£,r} then the (£, r)-cumulant of p.bP,a,p with respect to x is Aofl^ bf^Dlyjlx 

If, for each A^ G N, X^{N) = aL{X(N)) and X'^{N) = PR{X{N)), we claim the joint distribution of the 
random pair of matrices {X^{N), X^{N)) with respect to := ;^Tr o Ec{Xn) tends to pibP,a,p as N ^ oo. 
Indeed, by the moment-cumulant equations o and ([2]) it suffices to show for all x : {Ij ■ • ■ > t {.^, r} that 

Jim 

Notice that if Xf, Xr ■ {1; • ■ •: ''^1 {£, r} are defined by xe{k) = £ and Xr{k) = r for all k then 

Jim (xx(i)(Af),...,Xxi”i(Af)) = Jim K^iaX (N),..., aX (N)) = Xa^ and 

Jmr^ ..., = Jim^ k„(^X(1V), ...,f3X{N)) = Xp^ 

since 'd! m{L{X{N))'^) = dt m{R{X{N))'^) = dt m{X{N)'^) (the later computed traditionally). Therefore, if 
a = 0 or /3 = 0, the proof is complete. Thus we assume that a, /3 0 and x Y Xi-, Xr- 

Let s be the permutation such that 

X”^(W) = {s(l) < ■•■ < s(fc)} and x“^(W) = {s(* + 1) < ■ • ■ < sW}- 


Let X = X° s. Note that replacing x by X corresponds to moving all the right nodes down to be beneath the 
left ones, without changing their relative order. Furthermore s induces a natural isomorphism from BNC{x) 
to BNC{x) via tt i—>■ • tt such that lx) = t>‘BNc{s~^ ■ tt, 1^) as pbnc is completely determined 

by the lattice structure (see Definition 12.41) . Therefore, since X^{N) and X'^{N) commute, we obtain that 




= E 'I'N.(^''^'HlV),...,X>^(")(lV))AiBAc(^,lx) 

ir^BNCix) 

= V ^MAX\N),...,XYN),XYN),...,XYN))piBNc{Tr,l^). 

* ^ s. _ _✓ s._ _✓ 

TTGBNCix) 


|X“R{^})I copies 


lx“b{c})l copies 


Let x' ■ {^T ■ ■ ,rL} —>• {£,r} be defined by x^(n) = £ and x'{k) = x{k). Note that replacing x by x' 
corresponds to moving the bottom right node to the left side and thus induced a natural isomorphism from 
BNCix) lo BXCix')- Since X^{N)Im = —X^{N)Im, we obtain, as in Theorem 10.2.1], that 


= - V ■d!MAX\N),...,xYN),XYN),...,XYN),xYN))pBNc{T^Yx')- 

rv < * ^ ^ 


TveBNCix’) 


|X“T{^})I copies 


|X“^({i"})| —1 copies 
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By repeating the above arguments (i.e. commuting X^{N) past all the X‘^{N) in the above expression and 
then changing the last X‘^{N) to a X^{N)), we obtain that 

am 

Xa- = Aa'^ 

as desired. 

We will generalize a few more results from random matrix theory to the bi-free setting by using Theorem 
14.131 Of course these results are not exhaustive and there are many more applications and generalizations 
of random matrix results. 

For the first result, we will use the ideas of [101 Lecture 22] to demonstrate the asymptotic bi-freeness 
from Gaussian random pairs of matrices and “constant matrices”. By constant matrices, we mean elements 
of Mn{C) C Mn{^oo{^, fJ')), say {Dm{N)}meK', such that there exists a non-commutative probability 
space (Al, v^) and elements {dm}meK' C A such that the joint distributions of {DmiN)}meK' with respect 
to ;^Tr tend to the joint distribution of {dmImGif'- Often one takes Dm{N) to be diagonal matrices as such 
matrices can approximate the distribution of any measure on C. 

The following result is immediately obtain from the proof of nni Theorem 22.35] by running said argument 
with pairs of Gaussian random matrices independent from one another instead of just independent Gaussian 
random matrices - something that had little value to write down in the past. As the proof is modified simply 
by keeping track of cumulant terms from non-crossing pair partitions (e.g. see the proof of Theorem I4.10p . 
we omit the details. For those unhappy with this omission, the results of Sections [5] and [6] can be used to 
prove asymptotic bi-freeness from ‘constant diagonal matrices’ (see Remark 16.4|) . 

Theorem 4.16. For each N let {{Zk{N), Z'f.{N)y\k^K be pairs of N x N Gaussian random matrices 
such that {Zk(N), Z'f.{N)} is independent from {Zm{N), Z!^{N)} if k ^ m and the covariance matrix of 
{ZkiN), Z'f.(N)) is a positive, non-singular matrix Ck with real entries. For each iV S N let {Dm{Ny\m^K' Q 
M w(C) be as above. 

The joint distribution of 

{{ZkiN), Z'k{N))}keK U {Dm{N)U^K' 
with respect to ^Ty o tends to the joint distribution of 

b) {dm^mCiK^ 

where {{sk,Sk)}kGK is free from {dm}mGK', (sfc,Sfc) is independent from {sm,s'^) ifkf^m, and (sfe,sj,) is 
a pair of self-adjoint semicircular variables with covariance matrix Ck- Consequently, Theorem \4.1S\ implies 
{{L{Zk{N)), R{Z'f.{N)))}k^K are asymptotically bi-free central limit distributions and are asymptotically bi- 
free from {{L{Dm{N))}meK'ARiDm{N))}meK')- 

For our final application of Theorem 14.131 we will examine Haar bi-unitary random pairs of matrices. 
Recall from [21 Section 10] that a pair {Ut,,Ur) in a non-commutative probability space {A,^}) is said to be 
a Fdaar bi-unitary if the joint distribution of {Ui, Ur) with respect to q) is equal to the joint distribution 
of {U,U) where C/ is a Haar unitary in some non-commutative probability space. Gonsequently, if U{N) £ 
Af Ar(Loo(f2, /i)) is a Haar unitary random matrix, we will call {L{U{N)), R{U{N))) a Flaar bi-unitary random 
pair of matrices. 

Based on definition and uni Theorem 23.13], we immediately obtain the following. 

Theorem 4.17. For each N G N let {UkiN)}k^K be independent N x N Flaar unitary random matrices 
and let {Dm{N)}TneK' ^ A4n{<C) be such that for some *-non-commutative probability space {A, (p) and for 
some elements {dm}meK' ^ A the joint *-distributions of {Dm{N)}m£K' with respect to ;^Tr tend to the 
joint *-distribution of {dm}m^K' ■ Then the joint *-distributions of 

{{L{Uk{N)), R{UkiN)))}keK U {{L{DUN))U<,k'AR{DUN))U^k') 
with respect to -^Ty o Ec(^Xn) tends to the joint *- distribution of the bi-free pairs of faces 

{{U^^k: Ur,k){kCiK U {{L{drr){mCiKG {R{dm){mCiK') 
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as N ^ oo, where each {Ue^k,Ur^k) is a Haar bi-unitary. In particular, Haar bi-unitary random pair of 
matrices are asymptotically bi-free from constant matrices. 

Remark 4.18. Recall [U Theorem 10.1.3] shows that if (C, D) is a pair of algebras that is bi-free from a Haar 
bi-unitary {U^, Ur) in a non-commutative probability space {A, ip), then (C, D) and {U^CUf, UrDUf ) are bi- 
free pairs of faces. In particular, combining Theorem 10.1.3] with Theorem l4.171 by conjugating constant 
matrices by Haar bi-unitary random pairs of matrices, one can obtain many bi-free joint distributions. 


5. Bi-Matrix Models with Fock Space Entries 


In this section, a bi-matrix model involving left and right creation and annihilation operators on a Fock 
space will be examined. The results of this section generalize those of [121 Section 5] to the bi-free setting 
and provide bi-matrix models for all bi-free central limit distributions. 

Throughout this section let X = if (TL), namely the Fock space of a Hilbert space TL of sufficiently large 
size. If H is the vacuum vector of iF{'H) then px : d:” —>■ C is defined by pxi^ft (B f]) = X. Furthermore 
Po '■= Px '■ —>■ C is defined by po{T) = (TH, H) (all inner products linear in the first entry) thereby 

making {C{X),pq) a non-commutative probability space. 

Given an element h G H, let l{h) and r{h) denote the left and right creation operators by h respectively 
and let l*{h) = {l{h))* and r*{h) = {r{h))* be the corresponding annihilation operators. Then: 


l*{hl)l{h 2 ) = {h 2 M)Ix, (6) 

r*{hMh2) = {h2M)Ix. (7) 

[r(hi),r(h2)] = r(hi)r(h2) - r{h2)l*{hi) = (h2,hi)Po, (8) 

[r*{hi),l{h 2 )\=r*{hi)l{h 2 )-l{h 2 )r*{hi) = {h 2 ,hi)Pa, and (9) 

[l(/ii),r(h 2 )] = [r(hi),r*(h 2 )] = 0, (10) 


where Pn G C(X) is the projection onto the vacuum vector. Furthermore, for all T,S G C{X), 


piTPnS) = {TPnSn,n) = p{S){Tn,n) = piT)piS). 

Using notation and conventions from Construction IXTl we have the following which is the bi-free analogue 
of part of m Theorem 5.2]. 

Theorem 5.1. Given an index set K, an N G N, and an orthonormal set of vectors {h^j \ i,j ^ 
{1,..., N}, k G K} C P consider the N x N left and right matrices of C{Xm) defined for each k G K 
by 

Lk-.= ^L{[l{hl^)]), Ll-.= ^L{[P{hl,)]) Rk-.= ^R{[r{hl,)]), Rl := ^R {[r*{hi,)]) . 

Let^ = APvo Ec(Xn) ■ Then 

(1) L^Lk = Sk,mIxN, Rm^k = Sk,mIxM> [Lm, Rk] = [L^, R^] — 0 , [L^, Rk] = Sk,mPo, and [R^,Lk] = 

dk,mPo where Pq G £(Xj\i) is the linear map 


7"o(fe.i]) = -diag(Pa(Tr([^,,,])),..., Pn(Tr(K,,,]))). 


(2) = Po- 

(3) ^{TPoS) = $(r)$(S') for all T,S G C{Xn). 

(4) The joint distribution of {Lk, Rk, RkfkeK with respect to 4) is equal the joint distribution of 
{l{h^),l*{h^),r{h’^),r*{h^)}k^K with respect to po where {h'^fk&K UR is an orthonormal set. 


Proof. For m , notice 


L*mLk = 


N 




x—1 


jY ^k.,mL 


N 


x—1 
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as L is a homomorphism and 

1 /r ^ 

-y=i 


R*^Rk = 




■ N 

.y=i 


= Sk,mR{lN,C{X)) = 5k,mIxN 


as i? is a homomorphism on A4n{£{X)°'^)°^. Furthermore equation (ITUl) together with Remark I3.2I implies 
that [L^,Rk] = [L*^,RI]=0. 

Since 




N 


and since, by equation 


RlLn.[UA = -^Rl 


Ly=i 


N 




1 

TV 


N 




N 




N 

1 

TV 


— LraRkl^ij] + ^ 


.a:,y=l 

■ N 

N 

E ihyh^.y)Py^,y 

-x,y^l 

N 




^ ^ ^i,jPQ^X,3L 




we obtain that [L^^Rk] = Similar computations show [R^^Lk] = thereby completing part 

O- 

Part (m is a trivial computation using the fact that = Pq. To see part ([3]), write S{In^q) = 

Then 

1 ^ 

yy = j;^T.pyyy) 

y^l 

SO 

PoS{lN,n) = ^diag^Pn ^E^y.y^ ,---,Pn = ^'(5')/Ar,n- 

Hence 

$(TPo^) = pxj,{TPoS{lN,n)) = $(5W(T/A,n) = $(5)$(T). 

For part o, it must be shown that given any word W in {Lk, L^., Rk, Rl.}keK we have $(VF) = ipo{w) 
where w is the corresponding word in {T(/i^),r(h^),Note [T^ Theorem 5.2] completes 
the claim if no Rk nor PJ appear in W. Thus we proceed by an induction argument using parts o and ([3]). 

Suppose it has been demonstrated that <i>(kF) = <po{w) for all words W with at most n > 0 occurrences 
of {Rk, Rl}keK and all words W with n + 1 occurrences {Rk, Rl}keK and length at most to > 0 . Let W be 
a word with n + 1 occurrences of {Rk, Rk}k€K and length to + 1 . 

If bF = VR^ for some word V of length one less than the length of W and some g G {•,*}, then, since 
Rk^P = Lk^N,n (see Remark [33]), we obtain that 

$(bF) = ^{VRl) = ^{VLl). 

If ic = vr^ is the corresponding word in {l{h^),l*{h^),r{h^),r*{h ^)}then since we obtain 

that 

po{w) = po {vrD = ipo {vlD . 

Thus the inductive hypotheses then implies that $(IF) = ipq{w). 
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Otherwise we may decompose W = ViR^^L^V 2 for some 51,52 G {1, *} and some words Vi, V 2 with V 2 a 
word in {Lk, L^jk^K- We will now demonstrate a process to reduce the length of L^V 2 . Write w = virf}l!^V 2 - 
If 5 i = 52 or 5 i ^ 52 and k ^ m, then, by part o along with equations ( 0 , ®, and (nni), we obtain 

1 ^ 2 ) and (po{w) = ipo {viV^rfv 2 ) ■ 

Otherwise, if 51 = *, 52 = 1, and k = m, then 

^W) = $ {ViLS^RfV2) + HV 1 P 0 V 2 ) = $ {V,LS^RfV2) + $(b^i)$(1/2) 
by parts 0 and 0 whereas 

(po{w) = ipo {vil^rl^V2) + IPo{viPqV 2) = (fio ^^2) + ‘Po{vi)‘Po{v 2 ) 

by equation 0. As similar expressions hold for 51 = *, 52 = 1, and k = m (changing + to —), by 
the inductive hypotheses we obtain that $(bh) = (po{w) provided ^ {ViL^Rf^V 2 ) = (fo V 2 )■ By 

repeating this process, we eventually reduce the word W into one ending with R^^ thereby completing the 
inductive step by previous arguments. □ 

Remark 5.2. As all bi-free central limit distributions may be obtained using left and right creation and 
annihilation operators on a Fock space by m Theorem 7.4], Theorem IS.ll mav be used to construct bi-matrix 
analogues of all bi-free central limit distributions by taking suitable linear combinations. 

With Theorem 15.11 complete. we now demonstrate the bi-free analogue of the second half of [T^ Theorem 
5.2]. 

Theorem 5.3. With the notation in Theorem \5.1l 

{{Lk, Ll}keK, {Rk, Rk}keK) and (L(A4Ar(C)), i?(Af Ar(C)°P)) 
are bi-free with respect to $. 

Proof. To prove the claim, it suffices by Theorem l2.6l to show that the universal bi-free moment polynomials 
from equation 0 are satisfied; that is, we must demonstrate that 


$ 


/^X(l),e(l) . . . ^xW)An) 

\ 


E 

TreBNCix) 


X] h'BNc{'^,cr) 

o-eBNCix) 

7r<(T<e 






x(i).di) 

fcl > 


^X(n),e{n) 

■ • > 


( 11 ) 


for all X : {1,..., n} —>• {£, r} and e : {1,..., n} —^ {1, 2} where 

e{Lk,Lli;,,}k^K, e l{Mn{C)), g i?(Af^(C)°p). 

Suppose equality has been demonstrated in equation dill) for all y : {l,...,n'} —>■ {i,r} and e : 
{1,..., n'} —>■ {1, 2} with |x“^({r})| < m and for all y : {1,..., n'} —>■ {£, r} and e : {1,..., n'} —>■ {1, 2} 
with |y“^({r})| = m -I- 1 and n' < n. Note that the base case to = 0 follows from [TU Theorem 5.2]. Fix 
y : {1,.. .,n} {£,r} and e : {1 ,..., n} {1,2} with |y"^({r})| 

If y(n) = r, let y' : (1,..., n} —>■ {£,r} be defined by y'(n) = £ and x'iu) = xiv) for all y ^ n. Since 

RklN,n = LkiN.a, RIIn.q = LllN,n, and ])/Ar,n = L{[aij])lN,n, 


by replacing with Z^ obtain 






f 7X(l),e(l) r,: 

\'^ki 


X(n),e(n)^ 

kn 

X{n),t{n] 

kn 


) = $ and 

di) 


= 






x'(n).e(n) 


) 


where tt G BNC{x) is automatically an element of BNC{x') as moving the bottom node from the right side 
to the left side is a bijection from BNC{x) to BNC{xf). Since this bijection preserves the coefficient 


yBNc{'^,o-), 

a&BNCix) 

7r<tT<e 
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due to the relation between fiBNC and fj,NC (see comment after Definition 1^ . we obtain both sides of 
equation (ED are preserved under this operation. Consequently, equation ED holds in this case by the 
inductive hypothesis. 

Otherwise, we may select x € {!,...,n} such that x(cc) = r yet xiv) = ^ for all y > x. Let x' • 
{1,.. . ,n} —>■ {£,r} be defined by x'(x) = £, x'(a; + l) = r, and x'iu) = xiu) for all y ^ x,x + l. Similarly, let 
e' : {1,..., n} —>• {1,2} be defined by e'{x) = e{x + 1), e'(x + 1) = e(a;), and e'(y) = e(y) for all y x, x + 1, 
and let 


r/x'{x) {x) _ 7X(2:+l).e(a;+l) yx'{x+l) _ yx{x),t{x) 

K ~ ’ K+i ~ 


and 




■x'(y).e'(y) 


^x{v)Av) 

ky 


for all y x,x + 1. Note there is a bijection from BNC{x) to BNC{x') obtained by sending an element 
TT G BNC{x) to tt' G BNC{x') where tt' is obtained from tt by interchanging x and x + 1. 


If 


yx{x),e{x) yx{x+l),e{x+l) 


= 0 


then for all tt G BNC{x) 


$ 


('^xdl.dD ... and 






X(").e(n; 


= 


/ 7x'(l).e'(l) yx'{ri),e(n] 


where the second equality holds for all tt as either x and x+1 are in the same block of tt and the corresponding 
operators commute, or are in different blocks in which case the equality is trivial. Since 


y-BNc{T^,cr) = yBNcW,(k') 

creBNCix) cr'eBNCix) 

7r<cr<e 7r'<(T'<e' 


due to the connection between (Ibnc and fiNC, both sides of equation ED are preserved under this operation 
in this case. Otherwise if 


yX{x)M^) 7X(a;+l),e(2:+l) 


^0 


then we must be in the case 


or the case 


7Xix),e(x) 


= Rk and 


yx(x+l),e(x+l) _ J* 


z. 


X(x),e(x) 


= Rt. and Z, 


X{x+l),t{x+l) 


— Lk 


for some k G K. In the first case, we obtain that 


$ 


( zj !"' 


e(l) . . , 2;^(n),e{n 




$ 




(1) yX(=o-l),e(x-l) 

■ ■ ■ 


r ) (b f T yX(x+2),e(x+2) yX(n),e(n)\ 


by parts (HD and ([3]) of Theorem 15.II If x 9^',^ a; + 1, then we trivially obtain that 






e(l) yxin),e{n) 




However, if x a: + 1, then 


$ 

= $ 


tv' , • ■ • , j 

4- * ( 7 X(a^-l).e(a:-l) r \ ^ (j yx{x+2),t{x+2) 7 X(").e(n)\ 


where tti = 7r||-i ,j,} (remove all of x +1,..., n otherwise keeping the blocks the same) and 712 = 'x\{x+i,...,n}- 

Note the map taking tt G BNC{x) with x a: + 1 to (771,772) G BNC{x\{i,...,x}) x BNC{x\{x+i,...,n}) is a 
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bijection such that 






X] T-BNciA,(y) = 

TBNc{t^1,0-i) 


MBNc(t2, 0 - 2 ) 

aeBNC(x) 

aieBNC(x\ii ..}) 


cr2eBAC(x|{x + l,...,n}) 

7r<(7<e 



■7r2<o'2<e|{x+l,...,n} 


where cti = cr|{i,...,a;} and 02 = Indeed the above equality can immediately be obtained us¬ 

ing Definition 12.41 (which shows ^jlbnc is completely determined by the lattice structure) followed by the 
multiplicative properties of the Mobius function to obtain 

Ms Arc (t’’? O') = U7r2,fTi U(T2) = 


Since, the inductive hypothesis implies 


$ U: 


x(i).di) 

ki 


■ ■■z 


Xix-l),e{x- 
f^X — 1 


^ S> f T yX{x+‘2),e{x+2) 7 X("),e(n)\ 


E 

TrieBNC(x\{i,...,x}) 
Tr 2 &BNC(x\{x + l,...,n}) 







y^ MBArc(Ti, Ui) 


y^ MBArc(T2, 0 - 2 ) 

a-ieBNC{x\{i,...,x}) 


o- 2 GBAfC(x|{x+i,...,n}) 

7ri<CTi<e|{i,.,._x,} 




7X("),e(") 


/ 7X(1),«(1) ^x(x-l)Ax-i) T \ (j, (T ^x(x+2),e{x+2) 

we obtain that equation (fTTl) holds for the sequence ..., if and only if equation (ITTl) 

holds for the sequence (^Z^, ... Z^, provided 

^ ^xix+i)Ax+i) ^ 

The same holds by similar arguments when 


yX(x),e(x) _ p* 

'^kx — 


and 


By repeating these arguments of interchanging left and right operators (which preserve both sides of 
equation (fTT]) l we eventually reach the case that x(n) = r. As a previous case implies equation (fTTl) holds 
when x(n) = r, equation (fTTl) holds. Thus the result follows by the Principle of Mathematical Induction. □ 


6. Bi-Matrix Models with g-DEFORMED Fock Space Entries 

In this section, we will demonstrate the bi-free analogue of m Theorem 2.1] by determining the asymptotic 
distributions of left and right matrices containing creation and annihilation operators on g-deformed Fock 
spaces. Although this section is long and technical, it can be summarized via two points: 

(1) If {hk}k&K is an orthonormal set in a Hilbert space TL, then the two-faced families 

{{{l{hk),l*{hk)}Ar{hk),r*{hk)})}keK 

are bi-free and a bi-free central limit distribution by m Theorem 7.4]. Thus the joint distributions 
are completely determined by bi-non-crossing pair partitions where the *-terms precede the non-*- 
terms and the pairs must be of the same k G K colour. 

(2) The proof of [TTl Theorem 2.1] directly generalizes once the permutation is applied in the appro¬ 
priate manner. 

As a consequence of Theorem 16.11 Theorem 14.111 is automatically obtained via the g = 1 case and many 
implications from m immediately carry forward to the bi-free setting. 

We begin with definitions. Let TL he a (finite dimensional) Hilbert space. Recall, for g G [—1,1], the 
g-deformed Fock space is the Hilbert space 

nj®n 



Tq{Ti)=Cnq(B 
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equipped with the inner product 


N 


{gi 0***0 g-n-) 0***0 hm}q — ^n,m ^ ^ g ^ ^ {gkt 

<T^Sn k—1 

where Sn is the permutation group on {1 ,..., n} and 

mv(cr) = \{ii,j) €{l,...,ny \ i < j,a{i) > a{j)}\ 
is the number of inversions in cr. We define the vacuum state ipq : L{Tq{T-L)) —>■ C by 

For this section X will denote Fq{'H). 

For each h G H, there are left (right) creation and annihilation operators, denoted a{h) and a*{h) {b{h) 
and b*{h)) respectively, which are defined by 

a{h)Qq = h, a*{h)Qq = 0, 

b{h)Vlq=h, b*{h)nq=0, 

a{h){hi ® ■■■ ® hn) = h (g) hi (g) • • • (g) h„, 
b{h){hi (g) • • • (g) hn) = hi ® ® hn ® h, 


a*{h){hi 1^1 ■ ■ ■ <Si hn) = q’' ^{hk,h)hi igi ■ ■ ■ i^i hk <■ 

fe=i 
n 

b*{h){hi 0 • • • 0 h„) = ^ q^~^{hk, h)hi ®---®hk' 


hn, and 


fe=i 

where hk denotes hk is omitted in the tensor. Note {a{h))* = a*{h) and {b{h))* = b*{h). Furthermore, it is 
possible to verify the relations 

( 12 ) 

(13) 

(14) 


[a(hi),6(h2)] = 0, 

[a*(hi),6*(/i2)] = 0, 

[a*(hi),5(h2)] = (/i2,hi) , 

k n>0 


[ 6 *(hi),a(h 2 )] = (h 2 ,hi) I , 


(15) 


k n>0 


a*{hi)a{h 2 ) - qa{h 2 )a*{hi) = {h 2 ,hi)Ijr^(qi), and (16) 

b*{hl)b{h 2 )-qb{h 2 )b*{hi) = (h 2 ,hl)/^,(„), (17) 

where Pn G hl{X) is the projection onto the tensors of length n. 

Let IV G N and let K he a Hxed set. Suppose {h'^j \ z, j G {1,..., N}, k G K} is an orthonormal set in 
TL. For 6 G {^, r}, define 


zfiN) = ^Z^{[z^ihm) 


ZI’\N) = {[z\hl)]) 


1 


where, symbolically, = L, = R, = a, and z’’ = b. 

Theorem 6.1. The asymptotic distribution of 

{ZI{N),ZI’*{N),ZI’\N),ZI’**{N) I kGK,eG{i,r}] 
with respect to := o Ec(Xn) is equal to the distribution of 

z^’*, I k G K,9 G {£,r}'^ 
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with respect to ipo, where = I, = r, zl’' = z^{hk), zl’* = {z^)*{hk), zf,'* = z^{hl), = {z^)*{hl), 

and {hk,h\ \ k G K} is an orthonormal set. Furthermore, the difference in the distributions of a fixed word 
is O(^) as N ^ oo. 

Proof. This proof is obtained by modifying the proof of El Theorem 2.1]. As such, many details are omitted 
by referring to ra¬ 
ft suffices to show that ii x ■ ■ ■ ,n} ^ {I, r}, g : {1,... ,n} ^ {-, *, t, t*}, and k\,.. .k^ G K, then 

^Imi = ifio 


(18) 


Recall 


« £,- £,* £,t £,t* I 7 _ T^\ f r,- r,* r.t r.t* \ i _ \ 

.^fc \ k&K^,{z^ ,zf ,zf ,zf \kGK\'j 

is a bi-free central limit distribution by [M Theorem 7.4] with 

^ Afe' ) ~ 

unless k = k' and either w = * and w' = ■ ot w = t* and w' = t (in which case it is equal to 1). Consequently 
Fo = 0 unless there exists a tt G BNC 2 {x) such that 

(a) if X y then kx = ky, and 

(b) if X y with x < y, then either g{x) = * and g{y) = or g{x) = t* and g{y) = t. 

If such a TT exists, then cpo _ 7 i^y nioment-cumulant formulae of Section [5] 


By Lemma 13.71 and Remark 13.81 we obtain the left-hand-side of equation (fTBl) is 


1 


N 

Fq (^1 {js 


where 


i;fcl) ^2 {js-F2)’ds-\2) + F>^2)---Zn (j 

a{hl,) iig{m) = - 


TVf+ 1 ^ V''" V-^AAi) >^ 1 ) ^2 [Js-\2)^Js-\2) + F>^ 2 ) ■■■Zn („) + i 


;A:„)) (19) 


Zmii,j; k) = 


and 


Zm(,i,j; k) = 


o*(hfj ii g{m) = * 
aih^i) if gim) = t 
a*{hlj) iig{m) = t* 


Kkh) if gim) = ■ 

b*[h%) ifg(m) = * 

if gim) = t 

b*{hli) iig{m)=t* 


when x(to) = £ 


when xij^) = ^■ 


Notice that if h^j is replaced by zh^j for some z G C with jzj = 1, the distribution of any term in 


I'f j is replaced by zh^^ 

expression (HH) is not changed. Therefore, since h 1 —>■ a{h) and h 1 —^ h{h) are linear in h, we obtain that 

(^1 {ds^\ipJs-\l) + l:kl)---Zn (j.-b«)>A-b„) + i;fcn)) = 0 

unless there is a -tt G V 2 {n) such that 

(1) X y implies kx = ky, and 

( 2 ) X y and x y implies that the set {g{x),g{y)} is one {*, •}, {*, t}, {t*, •}, {t*, t}. 

Moreover, for such a partition tt, we have that if x ^,r y and x y then 

=js-\y)+i and = h-^y) if {gi.x),g{y)} = {*,-} or {U,t} 

7 - 1 / N = 7 - 1 / N and 7 - 1 / N , 1 = 7 - 1 / Nil otherwise. 

-'Sx m ■'^X (v) ■'^X A) + l ■'^X (y) + l 

Denote by 0 the set of all partitions satisfying © cLRci (|2p . Th 6 R GxprGssioR m equals 


( 20 ) 

( 21 ) 


1 


Xrf+i 


N 

E 


Fq [zi (j^-i(i), j^-i(i) + i; fcl) ■ • • Zn {is-\n)^ 3 s 


b") + l’ 


jk satisfy ]2QL J21i for some ttG© 
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The sum is zero if 0 is empty, which occurs when n is odd. Since tpo = 0 by above 

discussions when n is odd, we may assume n is even for the remainder of the proof. 

As in m, let C be a polygon with n edges and let 

D = {maps from the vertices on C to (1,..., N}}. 

Enumerate the vertices of C clockwise. Then the tuples {{ji, ■ ■ ■ ,jn) \ jk & , A^}} can be identified 

with elements of D in the trivial way where jk is the value of an element of D on the vertex. Orient the 
edges of C via 1—>-2—>••••—j-n—)-! and label the edge from fc to fc + 1 (where n + 1 = 1) with Sy^{k). Thus 
TT can be viewed as a partition on the set of edges of C labelled in this way. Consider the quotient graph 
C/tt, where two edges x and y with x ^ y and x y are “glued” together with orientation reversion if 
{gix), g{y)} = {*, •} or {t*, t} and with orientation preservation otherwise. Hence if s^{x) s^{y) then the 

vertices jx and jy+i pair up and the vertices jx+i and jy pair up if orientation is reversed and the vertices 
jx and jy pair up and the vertices jx+i and jy+i pair up if orientation is preserved. Note there is a clear 
injection i from 

Dtt := {maps from the vertices on C/n to {!,..., N}} 

into D such that (ji, ... ,jn) satisfies (1^ and (PT|) if and only if (ji, ■ ■ ■ ,jn) lies in the image of i : —>• D. 

Using this notation, expression m now becomes 


1 


ATf + l 


E 




for some ttG© 


(^i(. 




;fcl) (. 




; kn ^^ ■ 


As in the number of elements of Dt^ that are not injective is of order at most 0{N'i). Therefore, 
if D'^ is the subset of Dt^ consisting of injective functions then, modulo terms of 0(;^), the expression dT51) 
equals 

1 


iVf+i 


E 




for some ttG© 


(=■(■ 




;fci) (. 




; kn ^^ ■ 


If TT satisfies (P) and (P) and C/tt has at most ^ vertices, then |I? 7 r| = 0{N ^) so such tt may be ignored 
asymptotically in the sum. As C has n edges, C/tt has at most § + 1 vertices. By [TT] Lemmata 2.2, 2.3], 
such a C/tt has precisely f + 1 vertices if and only if tt non-crossing on (s;,^(l),..., s^{n)) (with the nodes 
in that order) and the edges of C were glued together only using orientation reversion. Hence tt G BNC{x) 
and we need only consider tt satisfying (EOD everywhere. 

Let 0 ' denote the subset of 0 consisting of tt G BNC{x) satisfying (1^ everywhere. It is possible to verify 
that \i{D'^) r\i{D/,)\ = 0{N^) whenever 7r,7r' G 0' are such that tt ^ tt' . Thus, up to 0{^), expression 

m equals 

E jvt+i E '^9 (^1 ■ 

Again, an empty sum is zero by definition. 

If TT G 0' is such that there exists x ^,r y with x < y and g{x) = ■ (so g{y) = *), then we claim that 


E (^1 (• 




;^i) 




; kn^ ^ 


= 0 . 


To see this, we recall that D/ consists of injective elements so the set equal to 

any set [js-\m)Cs^\m)+i} m^x,y. Since Zy {3s;i\yyjs^\y)+i^h) is either a*{K) or h*{h) for some 
h gH, we obtain using equations (IT^ . (fHl) . ([15]), (fTOl) . and dni that 
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for some to S N U {0}. Since Zy Vtq = 0, the claim follows. Similarly if tt S 0' is such 

that there exists x y with x < y and g(x) = t (so g{y) = t*), then 


(*■ (^- 


H ‘^9 

Thus we need only consider tt G 0' satisfying (0. 
Suppose TT G 0' satisfies ©■ We claim that 




kr,\]= 0. 


‘Pq 


(^1 ( j , 




; fcl) • • • (, 




; = 1 . 


Indeed, using just equations (HU), (USD, (HI, and m one can write 


( Js-l(l):Js-l(l) + l 


; fcl) • • • (, 




} kn^ 


as a product of 

i) 


i) 


where {h'l, ..., h'^} are unit vectors, g' : { 1 ,..., to} —>■ {•, +}, and there is a pair non-crossing partition 
tt' on {!,... ,to} such that if x y and x < y then g'{x) = *, g'{y) = •, and h'^ = h'y, and h'^l.h'y if 

X y, 

where {h},..., h!^ are unit vectors, g' : { 1 ,..., to} —>■ {•, +}, and there is a pair non-crossing partition 
tt' on { 1 ,... ,to} such that if x ~ 7 r y and x < y then g'{x) = *, g'{y) = •, and h'^ = and h'^l.h'y if 
X y, and 

hi) 

b*{h)a{h) and a*(h)b{h) 

where h is a unit vector. 

This is best seen via the following process on bi-non-crossing diagrams (thinking of equations (HU, (HU), (O, 
and m) where h^-Lky ii x ^ y. 


a*{hQ)> 


a*(/i7)T- 


a(/i 7 )^ 

a(/i6)f- 


b*ihs) 


-^b*{hs) 


^ 6 (^ 8 ) 


a(/i 5 ) *■ 
a*(hi) ' 
0 (^. 4 ) 


a 


-*b{hs) 

-^b*{h5) 

^h*{h2) 


a*{he) T 
a* (hr) 
a{h^)i 
a{he )i 


D 


-^b*{h^) 


,b*ih3) 


a{h^) f 

a* ( 114 ) 


a(/i 4 ) * 




Mh^) 

-*h{h2) 

-t&(/ii) 


n 


i&*(h 2 ) 

b*{h:i) 

b{h3) 

b{h2) 




-t6(hi) 


Since each of the four possible products when applied to Vtq produces fly, we obtain that 
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thereby completing the claim. Therefore 

whenever tt G 0 ' satisfies 

Hence, we have demonstrated that expression (fTOl) is equal to 

1 




;fcn)) = |i?; 


ATt + l 


E 


DL\. 


TreBNC2ix) 
satisfies lUi and 13 


The first part of this proof shows, either there are no such tt (in which case we get both sides of equation 
(USD equal to zero), or there is a unique tt (in which case the right-hand-side of equation (1181) was 1). In the 
case of a unique tt, since = 1 -|- O(-i) as D'^ is the set of injective maps on vertices, we 

obtain the right-hand-side of equation (TT^ is also 1 thereby completing the proof. □ 


As Theorem 16.11 parallels [TTl Theorem 2.1], we immediately obtain similar applications. Indeed asymp¬ 
totic bi-freeness from “constant block matrix algebras” holds. That is, for n a positive integer divisor of N, 
let 

M^{N)=Mn{C)®lK C Mn(C) <S> M^(C) C Mn(C). 

Theorem 6.2. Let 

{zf-(7V),zf’*(iV),4‘(iV),zf’‘*(iV) I keK,0e{£,r}j 
be as in Theorem \6.1\. and let n be a fixed positive integer. For multiples N of n, 


kGK 


is asymptotically bi-free from {L{Mn{N)), R{Mn{N))) with respect to Tjv := ;^Tr ° as N ^ oo. 


Proof. Each Z^'™{N) can be decomposed into an n x n block matrix with ^ x matrix entries. By 
Construction 13.ll Remark 13.51 and Theorem 16. 11 each ^ x ^ tends to a left or right creation or annihilation 
operator. As these left or right creation or annihilation operators appear in the correct entries to invoke 
Theorem 15.31 the result follows. □ 


Similarly, if A{N) C A^jv(C) denotes the subalgebra of diagonal matrices and {Dk{N)}keK C A{N) are 
diagonal matrices such that {Dk{N)}keK has a limit distribution as N ^ oo and sup^y^p^ ||Zlfe(lV)|| < oo for 
all k € K, then the following result holds as each Dk{N) can be approximated by elements of A{N) nT„(fV) 
for sufficiently large n and N. 


Theorem 6.3. Let 


{Zl--{N),ZI’*{N),ZI’\N),ZI’**{N) I kGK,9G{i,r}} 
be as in Theorem [Ql and let n be a fixed positive integer. Then 

{zi’-{N),zi’*iN),zi’\N),zi**iN)}^^^,{Zl'-iN),Zl'*iN),Zl’\N),Zl’**{N)} 


kGK 


is asymptotically bi-free from {{L{Dk{N))}k^K, {R{Dk{N))}k^K) with respect to at := o Ec,{Xf^) as 
N — >■ oo. 


Remark 6.4. We note that Theorem 16.II may be used to obtain many bi-free analogues of random matrix 
results as m Section 3] does. However, we will not formal state all such results. 

For one example, if g = 1 (the bosonic case) and {hk}keK U {^k}keK is an orthonormal set, then, as in 
m Section 3], a{hk) a*{hk) and b{hm) + b*(hm) are distributed with respect to ipq as Gaussian random 
variables of variance 1, commute with each other, and have a covariance matrix dependent on {hk,hm). 
Similarly a{hk)a* {h'^f) and b{hm) + b*{h'^^) are distributed with respect to ipq as centred complex Gaussian 
random variables of variance 1, commute with each other, and have a covariance matrix dependent on 
{hk,h'^) and {h'^,,hm)- Gonsequently, by using Z^’^N) Z^'*{N) for 6 G {£,r} and Theorem 16.11 which 

determines the limit distributions, several bi-free central limit distributions (contained in those where left 
operator commute with right operators in distribution) may be obtained. Furthermore, Theorem 14 .11 1 may 
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be obtain using this method and Theorems 16.21 and 16.31 imply such pairs of faces are asymptotically bi-free 
from “constant block matrix algebras” and “constant diagonal matrices”. These later results also hold in 
the fermonic case (g = —1). 

Similar results hold for “real Gaussian random matrices” by using the real and imaginary parts of Z^’’ (iV) + 
Z^’**{N) for 0 G {£,r} since Z^’'{N) + Z^’**{N) converges in distribution to a free circular element acting 
on either the left or right so that the real and imaginary parts are free semicircular elements. 

7 . Bi-Matrix Models producing Boolean and Monotone Central Limit Distributions 

In this section, we demonstrate the existence of left and right matrices of creation and annihilation 
operators on a Fock space that, when combined in certain patterns, have asymptotic distributions equal to 
Boolean independent, Boolean central limit distributions and monotonically independent, monotone central 
limit distributions. The main ideas behind these models arise from m, which demonstrates ways Boolean 
and monotone independence arise from bi-free pairs of faces. One motivation for these results is an alternate 
model for such distributions than those provided in [7]. 

Boolean Bi-Matrix Models. Recall, by an operator Z in a non-commutative probability space [A, (p) 
is a Boolean central limit distribution (with variance 1) if 

f 

Jv. 

where hb = \ (i^-i -I- Si) (the average of two point-masses; one at —1 and one at 1). Furthermore, operators 
Zi,..., Zn G A are said to be Boolean independent with respect to ip if 

m 

i=i 

whenever m G N, wi,..., Wm S N, and ki,..., km G {1,..., n} are such that kj ^ kj+i for all j. 

Let K be an arbitrary set and let {hj^k \ j > 1, fc G RT} be an orthonormal basis for a Hilbert space Ti.. If 
X = for each G N and k G K let 

LfJ L^J 

Tk{N) := l*{hj k) G E2j-i^2j{.N) -1- K^j,k) G E2j^2j+i{N). 

i=i i=i 

That is, Tfe is an upper triangular matrix with the sequence l*{hi^k), l{hi,k), G(li 2 ,fe), ^(^ 2 .fc),... on the sub¬ 
diagonal and zero elsewhere. Furthermore, let 

N 

S{N) ■.= J2Ej,j-i{N). 
t=i 

For each k G K let Lk{N) := L{Tk{N)) and R{N) := R{S{N)). 

Theorem 7.1. Using the above notation, for each k G K 

hm ^Tr ((Lfe(lV)R(iV))™)) = [ dp^B 

for all m G N. Furthermore, the operators {Lk{N)R(N)}k^K are asymptotic Boolean independent with 
respect to ;|Tr o Ec(Xn) ■ 

Proof. The idea of the proof is simply to compute the diagonal entries of 

L,(i)(lV)i?(iV)L,(2)(iV)i?(iV) • • • L,(„)(lV)i?(iV)J^,o 
for each n G N and e : {1,... ,n} ^ K. Indeed we claim the entry along the diagonal is 

f ^y,e(i)^y.«(2)^y+i,e(3)^y+i,e(4) \ix = 2y—\<N — n 

t ^y.e(l)^y+l_£(2)^i/+l,e(3)i ^y+ 2 ,e( 4 ) ■ ■ ■ H X = 2y < N — n 
y 0 otherwise 

This can be verified by appealing to Remark 13.21 and a direct computation of 

T,ii)iN)T,^2)iN)---T,iu)iN)S{N)-. 
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One then sees that 

^Tr {E^ (L,(i)(iV)i?(5(iV))L,(2)(iV)i?(5(iV)) • • • L,(„)(A^)i?(5(7V)))) 

_ 2 J max{0, if n is even and e(2m — 1) = e(2m) for all m 

iV \ 0 otherwise 

As limw-s-oo ^ niax{0, = 1, the result follows. □ 


Remark 7.2. Note one may instead use 


Ifl 1^1 

Tk{N) := ® E2j-i,2j{N) + ® E2j,2j+iiN) 

i=i i=i 

and the above computations will still work (i.e. we only need \K\ creation/annihilation operators). 


Monotone Bi-Matrix Models. For the monotone setting, as the main idea is taken from m where 
only monotonic independence of at most two algebras is obtained, we will only produce two asymptotically 
monotonically independent semicircular distributions. Recall from that the monotone central limit dis¬ 
tributions are also semicircular operators and two operators Z\ and Z 2 in a non-commutative probability 
space („4, (/j) are said to be monotonically independent with respect to ip if 

n+1 

ip (^z^^zi^zi^^ ■ ■ ■ = p Y[ V ■ 

i=i 


for all n G N U {0}, mi, m„+i G N U {0}, and kj,mj G N. 

Three types of operators will need to be considered. Again let X = T{'H). Fix A^ G N and let {/i„ | n G 
{0,1,..., A^}} be an orthonormal set in a Hilbert space %. Let s{hk) = l{hk) + l*{hk), let 


Ti{N) = L 




and 


Si{N) = R 


Te ] 

/ . Ejj+i \ . 


In addition, let 


T 2 {N) = L (diag(s(/ii),..., s(h„))). 

Let si,S 2 be two semicircular variables each of variance 1 that are monotonically independent with respect 
to a state ip. 


Theorem 7.3. With the above notation the joint distribution of Ti{N)Si{N) and T 2 {N) with respect to 
■Ew o Eik^Xn) asymptotically tend to the distribution of { 31 , 82 } with respect to ip. 

Proof. It suffices to show that 

hm lTr(F;£(;,^) (T2(Ar)™n7’i(^)^i(^))"^7^2(fV)’”^ • • • T 2 (iV)"^'‘(Ti(Ar).Si(iV))'="T 2 (iV)™'*+i)) 

N—^oo i \ 

— lb ( s'"! • • • s™” s''" 

- V (^S2 Si S2 Si S2 Si j 

where n,mi,mn+i G NU {0}, and kj,mj G N. 

Recall the definition of monotonic independence implies 

n+1 

ip (s^s^^sr • • • s+s^sr+^) = ip (sj^+-+''") n ^ (s?) • 

i=i 

We claim that 


T2 (A^)™i (Ti (Af)+(AT))''! r2(Ar)'"= • • • T2 (A^)"*" (Ti (Af)+(AT))''" T2 (A^)™"+i/iv.n 
is a diagonal matrix diag(2:i,..., zn ) where 

r 0 if X < kj 

“ \ s(+-^n)'"is(/io)''^s(/iy-i„_i)'"= • • • s(hy-ii)™"s(/io)''"s(+)™"+i if a; = 2/+ 
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Indeed the above can be verified since 


T 2 (iV)™diag(ei,..., = diag(s(hi)'"6, ■ • ■, s(/ia)™Ca) 

and 

(ri(Af)S'i(A^))diag(^i,..., = diag(0, s(/io)Ci, s(^o)6, • ■ •, s(/io)Ca-i)- 

for all diag(^i,..., ^a) G Xn- Since 


n+1 

ip {s{hy + r.r^s{h^f^ ■ ■ ■ s{hy+^r-s{h^f-s{hyr-^^) = p (s(ho)"^ + - + "") n ^ {s{hy + r. + l-ir^) 

t=l 

due to freeness, and since each s{hk) is a semicircular operator of variance 1, the result now follows as the 
zero diagonal terms are irrelevant as N ^ oo. □ 


Acknowledgements 

The author would like to thank Dimitri Shlyakhtenko for referring him to the non-commutative case and 
the papers [mill]- 


References 

[1] M. Bozjeko and R. Speicher, An example of a generalized Brownian motion, Comm. Math. Phys. 137 (1991), 519-531. 

[2] I. Charlesworth, B. Nelson, and P. Skoufranis, Combinatorics of Bi-Free Probability with Amalgamation, Comm. Math. 
Phys. 338 (2015), no. 2, 801-847. 

[3] _, On two-faced families of non-commutative random variables, to appear in Canad. J. Math. (2015), 26 pages. 

[4] K. Dykema, On certain free product factors via an extended matrix model, J. Fund. Anal. 112 (1993), no. 1, 31-60. 

[5] Y. Gu, H.-W. Huang, and J. Mingo, An analogue of the Levy-Hincin formula for bi-free infinitely divisible distributions 
(2015), 26 pp., available at arXiv: 1501.05369 

[6] R. Lenczewski, Matricially free random variables, J. Fund. Anal. 258 (2010), no. 12, 4075-4121. 

[7] _, Limit distributions of random matrices, Adv. Math. 263 (2014), 253-320. 

[8] M. Mastnak and A. Nica, Double-ended queues and joint moments of left-right canonical operators on full Fock space 
(2013), 28 pp., available at arXiv: 1312.0269 

[9] N. Muraki, Monotonic independence, monotonic central limit theorem and monotonic law of small numbers, Infin. Dimens. 
Anal. Quantum Probab. Relat. Top. 4 (2001), no. 1, 39-58. 

[10] A. Nica and R. Speicher, Lectures on the Combinatorics of Free Probability, London Mathematics Society Lecture Notes 
Series, vol. 335, Cambridge University Press, 2006. 

[11] D. Shlyakhtenko, Limit distributions of matrices with bosonic and fermionic entries. Free probability theory (Waterloo, 
ON, 1995) 12 (1995), 241-252. 

[12] _, Free Quasi-Free States, Pacific J. Math. 177 (1997), no. 2, 329-368. 

[13] P. Skoufranis, Independences and Partial R-Transforms in Bi-Free Probability, to appear in Ann. Inst. Henri Poincare 
Probab. Stat. (2015), 31 pages. 

[14] _, A Combinatorial Approach to Voiculescu’s Bi-Free Partial Transforms (2015), 19 pp., available at 

arXiv:1504.06005 

[15] R. Speicher and R. Woroudi, Boolean convolution, in Free probability theory. Fields Inst. Commun., vol. 12, American 
Mathematics Society, Providence, RI, pp. 267-279, 1997. 

[16] D. Voiculescu, Addition of certain non-commuting random variables, J. Fund. Anal. 66 (1986), no. 3, 323-346. 

[17] _, Multiplication of certain non-commuting random variables, J. Operator Theory 18 (1987), 223-235. 

[18] _, Limit laws for random matrices and free products, Invent. Math. 104 (1991), no. 1, 201-220. 

[19] _, Free probability for pairs of faces /, Comm. Math. Phys. 332 (2014), 955-980. 

[20] _, Free probability for pairs of faces II: 2-variable bi-free partial R-transform and systems with rank < 1 commutation 

(2013), 21 pp., available at arXiv: 1308.2035 

[21] _, Free probability for pairs of faces III: 2-variable bi-Free partial S-transform and T-transforms (2015), 16 pp., 

available at arXiv: 1504.03765 

Department of Mathematics, Texas A&M University, College Station, Texas, USA, 77843 
E-mail address: pskoufra@inath.tainu.edu 











